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Abstract 



We study the analytic structure of light-front wave functions (LFWFs) and its 
consequences for hadron form factors using an explicitly Lorentz-invariant formula- 
tion of the front form. The normal to the light front is specified by a general null 
vector The LFWFs with definite total angular momentum are eigenstates of a 
kinematic angular momentum operator and satisfy all Lorentz symmetries. They are 
analytic functions of the invariant mass squared of the constituents Mq = (J2 k^) 2 
and the light-cone momentum fractions Xi = ki-u/p-u multiplied by invariants con- 
structed from the spin matrices, polarization vectors, and These properties are 
illustrated using known nonperturbative eigensolutions of the Wick-Cutkosky model. 
We analyze the LFWFs introduced by Chung and Coester to describe static and 
low momentum properties of the nucleons. They correspond to the spin-locking of a 
quark with the spin of its parent nucleon, together with a positive-energy projection 
constraint. These extra constraints lead to anomalous dependence of form factors on 
Q rather than Q 2 . In contrast, the dependence of LFWFs on M5 2 implies that hadron 
form factors are analytic functions of Q 2 in agreement with dispersion theory and 
perturbative QCD. We show that a model incorporating the leading-twist perturba- 
tive QCD prediction is consistent with recent data for the ratio of proton Pauli and 
Dirac form factors. 



1 Introduction 

Light-front wave functions (LFWFs) are the interpolating functions connecting hadrons 
to their fundamental quark and gluon degrees of freedom in QCD. Many hadronic 
observables can be computed directly from these amplitudes. For example, matrix 
elements of local operators such as spacelike proton form factors, transition form fac- 
tors such as B — > ivir, and generalized parton distributions can be computed from 
the overlap integrals of the LFWFs. The determination of the hadron LFWFs from 
phenomenological constraints and from QCD itself is thus a central goal of hadron 
and nuclear physics. In principle, one can solve for the hadronic LFWFs directly 
from fundamental theory using nonperturbative methods, such as discretized light- 
front quantization, the transverse lattice, lattice gauge theory moments, or Bethe- 
Salpeter/Dyson-Schwinger techniques. Reviews of nonperturbative light-front meth- 
ods may be found in Refs. [TJ [2j [H| - 

One of the central issues in the analysis of fundamental hadron structure is the 
presence of nonzero orbital angular momentum in the bound-state wave functions. 
The evidence for a "spin crisis" in the Ellis- Jaffe sum rule signals a significant orbital 
contribution in the proton wave function [3J |5] . The Pauli form factor of nucleons 
is computed from the overlap of LFWFs differing by one unit of orbital angular 
momentum AL Z = ±1. Thus the fact that the anomalous magnetic moment of 
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the proton is not zero is an immediate signal for the presence of nonzero orbital 
angular momentum in the proton's LFWFs [Bj . It should be noted that orbital angular 
momentum is treated explicitly in light-front quantization; it includes the orbital 
contributions induced by relativistic effects, such as the spin-orbit effects normally 
associated with the conventional Dirac spinors. 

In this paper we shall show how orbital angular momentum is represented by light- 
front Fock state wave functions. A key tool will be the explicitly Lorentz-invariant 
formulation of the front form (see [1 for a review and references to original papers). 
The wave functions are defined at the light-front plane u-x = cr, for which the orienta- 
tion is determined by the null four-vector u. Although LFWFs depend on the choice 
of the light-front quantization direction lu, all observables such as matrix elements 
of local current operators, form factors, and cross sections are light-front invariants 
- they must be a;- independent. When the (^-independence is violated in approxi- 
mate calculations, one can still find the (^-independent form factors by separating 
them from the current operator matrix elements and omitting the non-physical, in- 
dependent contribution 7\. We shall show that the analytic form of LFWFs with 
nonzero orbital angular momentum is then constrained to a specific set of simple 
prefactors multiplying the scalar zero orbital angular momentum solutions. Know- 
ing the general form of the LFWFs can be important for determining the hadron 
eigenstates from QCD using variational or other methods. 

Our results for the analytic form of hadronic LFWFs, including orbital angular 
momentum, are consistent at large transverse momentum with the perturbative QCD 
counting rules of Ji, Ma, and Yuan |5J and the wavefunction constraints which fol- 
low the conformal properties of the AdS / CFT correspondence between gauge theory 
and string theory [TOJ IHj • 

We begin by noting that eigensolutions of the Bethe-Salpeter equation have spe- 
cific angular momentum as specified by the Pauli-Lubanski vector. The correspond- 
ing LFWF for an n-particle Fock state evaluated at equal light-front time a = uo-x 
can be obtained by integrating the Bethe-Salpeter solutions over the correspond- 
ing relative light-front energies. The resulting LFWFs ipi(xi,k±i) are functions of 
the light-cone momentum fractions Xi and the invariant mass squared of the con- 
stituents Mq = (J2i=i K) 2 = J27=i [ k±+m ] { and the light-cone momentum fractions 
Xi = k-u/p-u, each multiplying spin- vector and polarization tensor invariants which 
can involve The resulting LFWFs for bound states are eigenstates of a kine- 
matic angular momentum operator. Thus LFWFs satisfy all Lorentz symmetries of 
the front form JT2], including boost invariance, and they are proper eigenstates of 
angular momentum. 

There is now heightened interest in the analytic form of the nucleon form factors. 
Recent measurements of the proton form factors at Jefferson Laboratory ^Hl E| 
using the polarization transfer method, show a surprising result - the ratio 
G e{Q 2 ) / G m{Q 2 ) falls faster in momentum transfer Q 2 = —q 2 = —t than that found 
using the traditional Rosenbluth separation method. A possible source for this dis- 
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parity are the QED radiative corrections, since these are more likely to affect the 
Rosenbluth method f3 El EI- For example, the interference of one-photon and 
two-photon exchange amplitudes and the interference between proton and electron 
bremsstrahlung are present in the measured electron-proton cross section and can 
complicate the analysis of the energy and angular dependence required for the Rosen- 
bluth separation. 

If one translates the new polarization transfer results for Ge and Gm to the 
Pauli and Dirac form factors, the data appear to suggest the asymptotic behavior 
Q Fi(Q 2 ) / F 2 (Q 2 ) ~ const. In a recent paper, Miller and Frank ^H] have shown 
that the three-quark model for the proton LFWF constructed by Chung and Coester 
(CC) PH and extended by Schlumpf 2]Qj leads to QF l (Q 2 )/F 2 (Q 2 ) ~ const in the 
range of the JLab experiment, thus providing an apparent explanation of the JLab 
data. 

In dispersion theory form factors are analytic functions of q 2 , with a cut structure 
reflecting physical thresholds at timelike q 2 . This is also apparent from the analytic 
structure of Feynman amplitudes in perturbation theory. A functional dependence 
in Q = \/—q 2 does arise when there is a physical threshold at q 2 = 0, as in the 
case of gravitational 122] (or axial current) form factors due to the two-photon 
intermediate state; however, this would not be expected to occur for the vector current 
in QCD. 

Chung and Coester JH] introduced their ansatz for the form of baryon LFWFs 
in order to describe the static and low momentum transfer properties of the nucle- 
ons. The LFWFs in the CC model have the effect of spin-locking a quark with the 
spin of its parent nucleon, together with a positive-energy projection constraint. As 
we shall show, if one extends these forms to large transverse momentum, the extra 
constraints lead to an anomalous linear dependence of LFWFs in the invariant mass 
of the constituents and an anomalous dependence of form factors on Q rather than 
Q 2 . As we discuss in the conclusions, the lack of analyticity in Q 2 is related to the 
breakdown of the crossing properties incorporated in field theory. The CC constraint 
may provide a reasonable model for computing static properties of hadrons, but it is 
not applicable to large momentum transfer observables. 

We shall show that form factors computed from the overlap of LFWFs are analytic 
functions of q 2 due to their analytic dependence on the off-shell light-front energy and 
the general form of prefactors associated with nonzero orbital angular momentum. In 
particular, the general form of the LFWFs for baryons in QCD leads to a ratio of form 
factors F 2 (Q 2 ) / Fi(Q 2 ) which behaves asymptotically as an inverse power of Q 2 mod- 
ulo logarithms, in agreement with the PQCD analysis of Belitsky, Ji, and Yuan j'i.'jj . 
We also shall show that the form factor ratios obtained from the nonperturbative 
solutions to the Wick-Cutkosky model |2U have a similar behavior. 

It should be noted that the analytic form predicted by perturbative QCD is com- 
patible with the form factor ratio determined by the polarization transfer measure- 
ments. The detailed analysis of baryon form factors at large Q 2 based on perturba- 
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tive QCD predicts the asymptotic behavior Q 2 F 2 {Q 2 )/F 1 (Q 2 ) ~ log 2+8/(9/3) (Q 2 /A 2 ), 
where j3 = 11 — 2rif/3 [23J. This asymptotic logarithmic form can be generalized to 
include the correct Q 2 = limit and the cut at the two-pion threshold in the timelike 
region. Such a parameterization is 



F 2 /F 1 



l + {Q 2 /C 1 ) 2 \og b+2 {l + Q 2 /Aml) 
' l + (Q 2 /C 2 )Hog b (l + Q 2 /Aml) 



(1) 



where for simplicity we have ignored the small factor 8/9/5, as do Belitsky et al. For 
the large-Q 2 region of the available data, this already reduces to the asymptotic form 



F 2 /F 1 



K 



C|log 2 (Q 2 /4m 2 ; 



'C 2 



Q 2 



(2) 



Therefore, the values of C±, C 2 and b are not tightly constrained, except for the 
combination Cf/C 2 . A fit to the JLab data yields d = 0.79 GeV 2 , C 2 = 0.38 GeV 2 , 
and 6 = 5.1. Thus, as shown in Fig. Q one can fit the form factor ratio over the entire 
measured range with an analytic form compatible with the predicted perturbative 
QCD asymptotic behavior. 
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Figure 1: Perturbative QCD motivated fit to the Jefferson Laboratory polarization 
transfer data fSJ E]- The parameterization is given in Eq. ((TJ) of the text. The 
dashed line shows the predicted form for timelike q 2 = —Q 2 . 

The nominal 1/Q 2 power-law fall-off of F 2 (Q 2 ) / F\(Q 2 ) in perturbative QCD is 
a consequence of the underlying chiral structure of the vector interactions in QCD. 
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The factorized structure of hard QCD amplitudes predicts hadron helicity conserva- 
tion |2Zj at leading twist and thus the relative suppression of the Pauli form factor 
since it is a helicity-flip amplitude. These nominal power-law forms are also proper- 
ties of dimensional counting rules [231201130] for hard scattering amplitudes in QCD. 
The power-law suppression of i^Q 2 )/ Fi(Q 2 ) is not generally true for Yukawa the- 
ories with scalar gluons or in quark-diquark models of the nucleon based on scalar 
diquarks since the effective interactions violate chirality conservation. 

Iachello, Jackson, and Lande [3TJ IH21 EB] have introduced a model for the nucleon 
form factors based on dimensional counting and perturbative QCD at high momentum 
plus the analytic structure due to vector meson intermediate states. The model gives 
an excellent phenomenological description of the individual form factors and the form 
factor ratio measured using polarization transfer. 

Although the spacelike form factors of a stable hadron are real, the timelike form 
factors have a phase structure reflecting the final-state interactions of the outgoing 
hadrons. The analytic structure and phases of the form factors in the timelike regime 
are connected by dispersion theory to the spacelike regime. Each of the above models 
predicts a specific fall-off and phase structure of the form factors from s <-> t crossing 
to the timelike domain. As noted by Dubnickova, Dubnicka, and Rekalo, and by 
Rock [B3|, the phase of the form factor ratio Ge/Gm of spin-1/2 baryons in the 
timelike region can be determined from measurements of the polarization of one of 
the produced baryons in the exclusive process e _ e + — > BB, since the single-spin 
asymmetry normal to the scattering plane requires a nonzero phase difference between 
the Ge and Gm form factors. As demonstrated in Ref. [3H], measurements of the 
proton polarization in e + e~ —>■ pp will strongly discriminate between the analytic 
forms of the models which have been suggested to fit the proton Ge/Gm data in the 
spacelike region. The single-spin proton asymmetry is predicted to be large, of order 
of several tens of percent. 

The content of the remainder of the paper is as follows. The general construction 
of the wave functions appearing through the Fock decomposition of the state vector 
defined on the light-front uj-x = is explained in Sec. 12.11 The angular momentum 
properties of LFWFs are discussed in Sec. 12.21 For simplicity, we present the LFWFs 
for bound states of scalar fields and for a spin-1/2 system in Yukawa theory which 
can serve as a diquark-inspired model for a three-quark hadron. The consequences 
of the CC ansatz are explored in detail. Section E] recasts this discussion in terms 
of the standard light-front Fock-state expansion and the associated LF spinors. This 
allows contact to be made with perturbation theory. The construction of form factors 
and their asymptotic behavior are described in Sec. HJ The LFWF for the valence 
Fock state is obtained by integrating the Bethe-Salpeter solutions over the relative 
light-front energies. The expectation that LFWFs are functions of Mq, rather than 
Mq, is demonstrated in Sec. El The angular momentum properties of LFWFs are 
illustrated using the known nonperturbative eigensolutions of the Wick-Cutkosky 
model for nonzero angular momentum. SectionElpresents the implications for the CC 
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ansatz for a three-quark nucleoli state, where we show that the anomalous asymptotic 
behavior of QF 2 /Fx ~ const is a consequence of the extra constraints on the LFWFs 
imposed by the CC construction. SectionOcontains our conclusions, and an appendix 
collects some useful definitions and intermediate results. 

2 Light- front wave functions 
2.1 General construction 

The concept of a wave function of a hadron as a composite of relativistic quarks 
and gluons is naturally formulated in terms of the light-front Fock expansion at fixed 
light-front time, a = x-u ilj. The null four- vector u determines the orientation of the 
light-front plane; the freedom to choose u provides an explicitly covariant formulation 
of light-front quantization. For a stationary state we can consider a fixed light-front 
time and put a = 0. 

In this formulation, the eigensolution of a proton, projected on the eigenstates of 
the free Hamiltonian H^ D (g = 0), has the expansion 

M) = / ? / , A 1 ,...,A„(^'---' fc «'P' WT ) 

n>3 

^fcj-p-wr 2Mdr II (2 ^32 £fc 10) • ( 3 ) 



Here a) is the usual creation operator and = ymf + k 2 . All the four- momenta 
are on the corresponding mass shells: fc| = mf , p 2 = M 2 , (lot) 2 = 0. The set of light- 
front Fock state wave functions x (ki, . . . , k n ,p, lot) represents the ensemble of 
quark and gluon states possible when the proton is intercepted at the light-front. The 
Xj label the light-front spin projections of the quarks and gluons along the light-front 
quantization direction. The scalar variable r controls the off-shell continuation of the 
wave function. From the point of view of kinematics, the four-momentum lot can 
be considered on equal ground with the particle four-momenta k\, . . . , k n ,p. Being 
expressed through them (by squaring the equality J2] kj = p + tor), r reads 

M 2 - M 2 
2uj-p 

where Mq = kj) 2 . The difference M — M ~ r between the effective mass M Q of 
constituents and their bound state mass M is just a measure of the off-energy- shell 
effect. 

If the wave functions describe a system with spin composed of constituents with 
spin, they are represented through scalar functions multiplied by invariants con- 
structed from the spin matrices, polarization vectors, and In general, the scalar 
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functions depend on a set of scalar products of the four- moment a ki, . . . , k n ,p, lut with 
each other. One should choose a set of independent scalar products. A convenient 
way to choose these variables is the following. We define 

ui-p 

where S" =1 Xj = 1 and Y^=\ Rj — an d represent the spatial part of the four- 
vector Rj as Rj = R\\j + k±j, where R\\j is parallel to Co and k±j (with Yjj k±j = 0) is 
orthogonal to Co. Since Rj-uo = RqjUJq — R\u-oj = by definition of Rj, it follows that 
Roj = \R\\j\, and, hence, k\ - = —R 2 and k±rk±j = —R^Rj are expressed though the 
squares and the scalar products of the four- vectors Rj. Hence, they are the Lorentz 
and rotation invariants. Therefore, the scalar functions should depend on Xj, k\j and 
k±i-k±j. In terms of these variables, the integral in Eq. Q is transformed as 



„ / n \ n 

y...*W \Y,kj-p-u>T\ 2{u-p)dr n 



i (2tt) 3 2£ A 



n lit \ I \ 

J. ..25 E^-iM (2) e^- n 



n j2 



d }\i i jdoCj 



I \ / i 167T 3 Xj 

J / \ 3 ) »=1 * 

In this way we find that ip reduces to tjj^ Xn (xj, k±j). These wave functions are 
independent of the hadron's momentum p + = p-u and p±, reflecting the kinematical 
boost invariance of the front form. All observables must be invariant under variation 
of as we shall show, this generalized rotational invariance provides an elegant 
representation of angular momentum on the light front. 

A scalar bound state of two spinless particles with p 2 = (k a + k^ 2 = M 2 can 
be described using the covariant Bethe-Salpeter function $(xi,X2,p), which in mo- 
mentum space corresponds to &Bs(k a , k^) = &(q,p), where q = {k a — kb)/2. For 
J = the Bethe-Salpeter function is a scalar function of the Feynman virtualities 
with k% + k% =p». 

The explicitly covariant version of light-front dynamics lj turns into the standard 
one at the particular value of u = (1, 0, 0,-1). The conjugate direction is defined as 
(, where (-u = 2 and ( = (1, 0, 0, 1) in the standard light-front frame. In general, the 
two-body scalar LFWF can only be a function of the corresponding off-shell light-front 
energy 

M 2 -M 2 = (£ h-C - p<){uJ-p) = 2(u-p)r (4) 

i 

and the momentum fractions Xi = ki-uj/p-uj. 

One can identify the corresponding two-parton LFWF by calculating the Fourier 
transform of X2,p) with the arguments x±, X2 constrained to a fixed light-front 
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time cr = lu-x, i.e., the transform of the function $(xi,X2,p)5(u-Xi)5(u-X2)- In mo- 
mentum space it corresponds to integration over the minus component of the relative 
momentum k, or in the covariant form pQ 

^wwr }(HhpW . ( 5) 

n(u-p) J-oo 

The resulting LFWF describes a spin-zero bound state of on-shell partons with k 2 a = 
m l k b = m l 

The main dynamical dependence of a LFWF thus involves Mq = (J2ki) 2 , the 
invariant mass squared of the partons. Since the equation of motion for the LFWF 
involves the off-shell LF energy p-Q, the natural analytic dependence of (p is in the 
variable Mq, not M . 

It is useful to examine the properties of the LFWF in the constituent rest frame 
E?7] . We will distinguish the four vector = X) K and the bound-state four momen- 
tum p^ . Since there is no conservation law for the minus components of momenta, 
we have (-p 7^ and the constituent rest frame (V = 0) and the rest frame (p = 0) 
differ from each other. In the constituent rest frame where V = J2i ki = 0, we can 
identify 

Mo = p = E^ = Ev / ^ + < 

i i 

It is also convenient to make an identification of the front-form and the usual instant- 
form wave functions in this frame. 

In general, LFWFs are eigenstates of the LF angular momentum operator pQ 

J = — i[k x d/dk] — i[n x d/dn] + -a, (6) 

2 

where n is the spatial component of u in the constituent rest frame (V = 0). Although 
this form is written specifically in the constituent rest frame, it can be generalized to 
an arbitrary frame by a Lorentz boost. 

Normally the generators of angular rotations in the LF formalism contain inter- 
actions, as in the Pauli-Lubanski formulation; however, the LF angular momentum 
operator can also be represented in the kinematical form © without interactions. 
The key term is the generator of rotations of the LF plane — i[n x d/dn] which re- 
places the interaction term; it appears only in the explicitly covariant formulation, 
where the dependence on n is present. Details can be found in JT]. The application of 
the LF angular momentum operator © to the scalar function <p(Mq,x) verifies that 
it describes a J = state. 
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2.2 The Yukawa model 

The general form of the LFWF of a spin-1/2 system composed of spin-half and spin- 
zero constituents is given by 

ip(h,p) = u{h) [vi + — <^ 2 ] u(p), (7) 

where u = (uq,uj) is the four-vector determining the orientation of the light-front 
plane lo-x = 0. Here Cj = 0^7^, and u(k\) is the conventional Dirac spinor of the 
spin-half constituent, and u(p) is the spinor of the bound state. 

The general wave function (jJJ) in Yukawa theory is determined by two scalar 
components tp\ and </?2> each a function of Mq and x. In contrast, the wave function 
corresponding to the CC ansatz used in [THj is equivalent* to the form 

V> = c k c p u(ki)A + u(p)fi, (8) 

where P = £4 + &2 is the sum of the constituent momenta, M = v?^, = 



l/y/m + e k , c p = 1/yM + Ej,, e k = ym 2 + k 2 , and e p = \/M 2 + p' 2 (for simplic- 
ity we assume that constituents have equal masses m a = = m). The matrix A + 
in (jHJ is the projection operator 

A + = (9) 

+ 2M V ; 

which in the constituent rest frame becomes 

a + =(;°). do) 

Hence, in the constituent rest frame we obtain 

c k u(k 1 )A + = x\(l, 0), c p A + u(p) = ( J xn, (11) 

with Xi an d Xtv being two-component spinors for constituent 1 and the nucleon, 
respectively. From (jHJ) we reproduce the analog of the three-quark CC ansatz used in 
Ref. [IB], which is 

V£ = X t<T1 X CT /i = 5^h. (12) 
This should be contrasted with the general form (in the constituent rest frame) 

< = X t<T1 (/1 + £[n x / 2 ) X% (13) 



t Equivalence has been established [SSj in the sense that equivalent results are obtained for form 
factors. 
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where k = k\ = — k 2 and n = oo/\uj\ in this frame. The additional f 2 term represents 
a separate dynamical contribution with orbital angular momentum, to be contrasted 
with the purely kinematical contributions of orbital angular momentum from Melosh 
rotations. The form (fTB"J) can be found by substituting into Eq. the explicit form 
of the Bjorken-Drell (BD) spinors, as given in the constituent frame by Eq. (|59j) of 
the appendix. The two functions f\ and f 2 are then determined in terms of <pi and 



y?2 as 



fi = ^(M + m)[m(M + m) + M (xM + (l-x)M )]<p 1 

+ ^M(M + m) (xM + m)ip 2 , 
Mo 

h = ^ M ( M o + M )^-^(M 2 -M 2 V 1; (14) 
with x = xi = 1/2 — k-n/Mo- The inverse relations read 



(M + M)f x ~ tSfHM) + M) (xM + m) f 2 



CfcCp 



y 4M M V JJ 

+ aTLi ( M " + M ) t m ( M ° + m ) + M o( xM + - a; ) M o)] ^ • (15) 
4M Mfc 

Note that the wave functions (pi and ip 2 are Lorentz scalars and analytic functions of 
M5 2 and x. The derived amplitudes /i and f 2 contain kinematic factors linear in Mq 
which arise from the reduction of the covariant form to the Pauli spinor form in the 
constituent rest frame. The application of the LF angular momentum operator © to 
the LFWFs of the Yukawa theory in Eq. (|13|) verifies that these wave functions are 
in fact states with J = 1/2, J z = ±1/2. 

In the non-relativistic limit, where M = 2m, Mo = 2m, and x = 1/2, we get 

Vl = 2~7^ fl ~ V2 = W2 h ' (16) 

The component f 2 is of relativistic origin. In the non-relativistic limit c — > oo (here 
c is the speed of light), the LF plane t + z/c = turns into t = 0. Any n dependence 
in (fTHjl should then disappear, and this happens if f 2 becomes negligible. 

We see that in order to obtain the constrained CC form ()12|) from the general 
form (J7J), one should eliminate, according to (JTTJ), the second components of spinors 
and neglect the second component f 2 in (|13|). This is natural in a non-relativistic ap- 
proximation, when both the second components of spinors and the dependence on the 
light-front orientation indeed disappear. One can use this form of the wave function 
to estimate the influence of relativistic effects on the static nucleon properties j!9j . 
However, there is no compelling reason to use the form (|T2*|) in the asymptotic rela- 
tivistic domain. 
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The form ([12)1. determined by the single component fi, implies a relation between 
(fi and if2 in (JTj). Setting / 2 to zero in ([To]) , we obtain 



and, hence, 



= lt (M ° + M)A ' 



M -M 



Since M is large for large momenta, we see that in the asymptotic regime the com- 
ponent if2 dominates: <y9 2 ~ Moipi/(2M). 

Thus the ansatz (|12|) is equivalent to an assumption that the component if 2 dom- 
inates in (J7J). We shall show below that this predominance of if 2 over tpi, generated 
by the wave function (|T2*)I. results in the QF 2 /F 1 — > const asymptotic behavior. 



3 Light-front spinors 

The discussion of the previous section can also be given in terms of a LF Fock-state 
expansion and the associated LF spinors. The LFWF of a hadron with spin projection 
J z = ±| is represented by the function ip^ Xn ( x h ^-u)> where 

h = (*+ K, h±) = ( [ x l P + , k \^ p f\ h?j (19) 

specifies the momentum of each constituent and Aj specifies its light-front helicity 
in the z direction. The light-front fractions Xi = kf / P + are positive and satisfy 

= 1- We note that Mq = J27=i k±l ^ mi = h) 2 is Lorentz invariant, and the 
scalar part of the LFWF is a function of only Xi and M5 2 . 

For a spin-1/2 state with two constituents in Yukawa theory, we write xp^ z (x, k±) = 
ipj^ixi, k±i), where A = Ai is the helicity of the fermion, x = Xi, and k± = ki±. (We 
use the subscript 1 for the fermion and 2 for the scalar constituent.) The two-particle 
Fock state with total momentum (P + , P±) and spin J z is then given by 

P + ,P ± = U J Z )= f dX f 2k± Y,ti'(x,k A .)\xP + ,k ± ,\). (20) 
J l6ir s Jx(l — x) a 

The Fock-state ket on the right is defined by 

\xP + , k x , A) = \kf = xP + , k+ = (1 - x)P + ; k 1± = k ± , k 2± = -k ± ; A x = A) (21) 
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and normalized by 



k' iX , \'\k+, k i±} A) = J] lQv 3 ktS(k'+ - kt)5(k' i± - ku)5 x , tX . (22) 



»=i 



The four functions ^±1/2 anc ^ ^±1/2 provide a representation of the LFWFs for 
Jz — Tj J. POj- The associated spinors are the LF spinors u LF (101 E], which in the 
hadron rest frame (p = 0) are given by 

u LFai (k\) = 1 -.x Wl (m + k 10 + a z (a-h) - (k 10 - m)a z - a-h) , 

/2(fc 10 + ku) V ' 



The general form (JJJ) has a simple form in the light-front spinor representation. A 
straightforward calculation gives* 



VI - x^iii^, k±) = —<pi 



_ J 

2 X 



Vl-x^\(x,k±) = —<pi, 

+ 2 X 

y/l^x ^(x,^) = (M+™^cp 1 + 2Mcp 2 . (24) 
Note that in the perturbative Yukawa model§ one obtains 

^ 1= M^Mf ^ = °- (25) 

In this way we reproduce the wave functions (44) and (46) of [HH]- The general solution 

in (j24)L which does not require any assumptions, differs from the perturbative solution 

only by the contribution 2M(f2 in the components ■?/>!_ 1 and 

~'~ 2 2 



*The wave function 10 and the one given here are related as follows: ijj(m 00 = 
y/ 'x(l — x)^(here). In this expression, k x and k y are equivalent to k 1 and k 2 , and J z and L z 
denote J 3 and L 3 . 

Hil the case of perturbative models, a single-particle wave function V>£ = VZ5 2 {k ± )S(l- x)S JzXl 
is present, where the normalization constant Z ensures unit probability. The perturbative Yukawa 
model wave functions can be formally differentiated with respect to the boson mass in order to 
simulate the fall-off of the wave function of a composite hadron and eliminate the single-particle 
Fock component. 
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In terms of the LF spinor representation, the CC ansatz (fT2*|) implies particu- 
lar forms for the components of Substituting the expression (fT%|) for <p 2 into 
Eqs. (J24J), we find 



Vl -x^ +1 _(x,k ± ) = (m + 

r, , T / f \ ( + + iky) 



m 



x 



X 



R \i i T \ \~\~k x iky) 

VI - x ipl 1 (x, k±) = y —ifi, 

+ 2 X 
i _» / Tfl\ 

VT^Vlifofci) = { M o + -)Vi- (26) 

We see that the CC ansatz in (J12)) is equivalent to the replacement M — > M in the 
perturbative LF components given in [39j . The anomalous dependence on M is the 
source of the discrepancy of the ansatz with the asymptotic behavior of the compo- 
nents with different angular momentum projections found in jH] from a perturbative 
model based on the iteration of the one-gluon exchange kernel. The perturbative QCD 
counting rules for hard scattering exclusive amplitudes [2Hl E3 IHDj and the fall-off of 
hadronic LFWFs jH] at high transverse momentum can also be derived fOl IHj with- 
out the use of perturbation theory using the conformal properties of the AdS/CFT 
correspondence between gauge theory and string theory [TT] . 

The component ip_! corresponds to L z — 1, and the component ip i corresponds 

2 + 2 

to L z = 0. The PQCD analysis based on the exchange of a gluon gives [S] 

i>{L z = 1)/^{L Z = 0) = Vli/Vii ~ ~ l/k ± , (27) 

2 '2 

which implies Q % F 2 jF\ — const jH]. We confirm this in the next section. On the other 
hand, (fTHj) or (|2l)|) . corresponding to the CC ansatz, gives 

i>{L z = l)/i/>(L z = 0) = ^liMi ~ fc±PiM ~ const. (28) 

2 '2 

In turn, the ratio (J2HJ) results in the asymptotic behavior QF 2 / ' F\ = const, as we will 
see in the next section. 



4 Form factors and the Mq dependence of light- 
front wave functions 

4.1 Dirac and Pauli form factors 

In the case of a spin-| composite system, the Dirac and Pauli form factors F^q 2 ) and 
F 2 (q 2 ) are defined by 

(P'|J M (0)|P> =u{P') [F^Y + F^q^a^uiP), (29) 
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where g M = (P' — P) M , u(P) is the bound-state spinor, and M is the mass of the 
composite system. In the light-front formalism it is convenient to identify the Dirac 
and Pauli form factors from the helicity-conserving and helicity-flip vector current 
matrix elements of the J + current component jEj 



P + q,1 



P + q 
J+(0) 



J+(0) 



2P^ 



P,n=Fi(q 2 



2P^ 



p,l) = -(q 1 -«i 



2,F 2 {q 2 ) 



2M 



(30) 
(31) 



We use the standard light-front frame (q^ = q° ± g 3 ) where 

J2 \ 



q = (q ,q ,q± 
p = 



o 



-q 



' p j 



q± 



p + ,p-,p ± 



{ , M 2 - 



(32) 



and q 2 = —2P-q = 
photon. 

Using Eqs. (jSDJ and (jUJ) we have 



gl 2 is the square of the four-momentum transferred by the 



Fi(q 2 ) 



d 2 k±dx 
16?r 3 



where 



k r x — k± + (1 — x)(fj 
From Eqs. (|3*T|) and (j27Ijl . we have 



(33) 
(34) 



F 2 (q 2 ) 



-2M 



d 2 k\ dx r 



V>L*(zA)VdiM±) +^J_{x,k' ± )^K{x,k ± )]. (35) 



(fe - %) •> 1671-3 
The individual wave functions are given by (j2U); substitution yields 



FM 2 



d 2 k i dx 



167T 3 7 x 2 (l — x) 



(36) 



1 



x(2mM + xM 2 ) + m 2 +k 2 ± - -(1 - x) 2 Q 2 
+ 2Mx{m + iM)(^> 2 + if'^) + AM 2 x 2 ip' 2 ip 2 ) 



F2(q 2 ) 



M 



d 2 k±dx 
x 2 (l — x) 



(37) 



x 1(1 - x) (m+ xM) ip'^i - 2Mx kj ^' L ((pi<p' 2 - (p2(p[] 

+Mx(l - x)(<fi<f' 2 + <f2<f'l)} , 
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where the kj_ variable has been shifted by — |(1 — x)q±, so that 

<Pi,2 = <pi,2{x, k_L ~ ^(1 - x)q±), <f' lj2 = <fi,2{x,kx + g(l - x)q±). (38) 

The required parton invariant masses M and M' are 

(k± - 1(1 - x)gl) 2 + m 2 (k_ L + Ul-x)q_ L ) 2 + m 2 

M 2 = A 21 J , M'l — ± 2 \ LJ . (39) 

x(l — x) x(l — X) 

Being expressed in terms of the invariant masses, the form factor integrands contain 
squares of masses: M5 2 and M'\. 

One can see that if (p 2 is smaller than <pi or of the order of (f\ for large k\, 
the leading terms in F\ are k\ and — 1(1 — a;) 2 <5 2 . An analytical calculation of the 
asymptotic behavior of the form factors with the power-law wave function (pi = 
N/(M$ + f3 2 ) n shows that these two leading contributions cancel each other. Thus 
the leading term becomes ~ logQ 2 instead of ~ Q 2 and the ratio F 2 /Fi becomes 
~ 1/ \og(Q 2 /m 2 ) . The same result is found for pseudoscalar coupling, i.e., with the 
wave function which is obtained by inserting in (J7J) the matrix 75. As explained in 
the introduction, the cancellation of the leading term in the Dirac form factor in the 
scalar gluon models is related to the violation of chirality conservation. 

We summarize in Table Q how the asymptotic behavior of the form factor ratio 
depends on the asymptotic properties of the LFWFs in the Yukawa model. Since the 
scalar part of the LFWF is a function of Mq and Xj, and ip(L z = ±1) contains the 
k x ± iky prefactor, ip(L z = 1)/i/j(L z = 0) can only be an odd power of k±. The third 
and fourth columns of Table 1 correspond to the cases of and ()27|). respectively. 



feg g k ± const I I 
F2 _L _L _L 1 J_ 

F! s W Q W s 



Table 1: The dependence of the asymptotic form of the form factor ratio on the 
asymptotic behavior of the LF components ip(L z = l)/ip(L z = 0) = ^li/Vli in the 

Yukawa model. 

The cancellation of the leading power-law contribution is specific to the scalar and 
pseudoscalar diquark models. In the case of a spin- 1/2 system comprising a spin-1/2 
quark and a spin-1 diquark, the hadron wave function is determined in general by six 
independent components. To see the effect coming from the spin-1 diquark, consider 
a wave function in the one-component form 

^a 1 ,x(k 1 ,k 2 ,p) = e^ x \k 2 )u ai (k 1 )^ u u(p)(p 1 . (40) 
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Here e^\k 2 ) is the spin-1 polarization vector. The sum over polarizations results in 
the propagator (<?„„/ — k 2u k 2u / / fi 2 ) in the form factor calculation. In this model we 
find for the form factor ratio 

F2 log(Q 2 /m 2 ) 

f 1 ~ q 2 

which is close to the fit (|2J), differing only by a factor log(Q 2 /m 2 ). The same per- 
turbative behavior (up to a coefficient) also occurs in the asymptotic behavior of the 
electron form factor in QED. If a scalar or pseudoscalar coupling is also present, the 
vector contribution will dominate the asymptotic behavior. We emphasize that for 
all three couplings, both form factors decrease as an integral power of Q 2 (modulo 
powers of \og(Q 2 /m 2 )), not as a power of Q. As shown in Fig. a fit to the form 
factor ratio based on powers of Q 2 and powers of log Q 2 describes the experimental 
data well. 



4.2 Consequences of the CC constraint 

The CC ansatz f 2 = is equivalent to a quark-scalar-diquark model, but with the 
additional condition (|18|) . As we have seen this introduces anomalous terms in the 
LFWFs which are linear in M . Substituting ip 2 from ([T%]) into Eqs. (JHtjj). (|3*T|). we 
obtain: 



1 



327T 3 J x 2 (l -x) 
(1 - x) 2 Q 2 l <p' m 



x(l - x)(M 2 + M'l) + 2x 2 M Q M' + 2xm(M + M' Q ) 



M 



d 2 k 1 dx 



167T 3 J x 2 (l — x) 



[2(1 - x)m + x(l - x)(M + Mq) 



x 



q2~(Mq — Mq) (M + Mq 



(41) 



(42) 



The same result is of course obtained by direct calculation with the wave function 



Because of the terms with Mq, 



M' , and MqMq, the form factor F\ contains 
the second power Q 2 relative to the term 2(1 — x)m in F 2 , which does not contain 
Mq or Mq. This occurs independently of whether the cancellation between k\ and 
— — x) 2 Q 2 takes place or not. If it does not take place, we get an extra contribution 
to the Q 2 power. This only changes the coefficient of Q 2 . Similarly, because of the 
term with (M + Mq), F 2 contains the first power of Q. This extra factor of M (and, 
hence, Q) results in the nominal asymptotic behavior 



QF 2 
Fx 



const. 



(43) 



17 



We emphasize that this behavior follows from the CC ansatz (fT2"]) (or f 2 = in ([13))) . 
which in turn, is equivalent to the relation ([T%]) . After the constraint is imposed, the 
asymptotic behavior fi3|) follows without further dynamical assumptions. 

The coefficient functions tpi are in general functions of Mq since in the equations of 
motion only the quantity Mq appears. In the next section this is shown explicitly in 
the case of the Wick-Cutkosky model. Explicit calculation of the form factors, with 
ipi = N/(Mq + /? 2 ) 3 5 , as assumed in Ref. [15] . and with ip 2 fixed by Eq. (JTHJ), qual- 
itatively replicates the nonconstant, nonintegral asymptotic behavior found in (18j . 
Figure E] shows the result for QF 2 /Fi, which is approximately constant for a small 
range of intermediate Q 2 but is well fit by Q~ 025 for large Q 2 . Similar results are 
obtained with an exponential wave function. 



1.2 




Figure 2: Asymptotic behavior of the form factor ratio QF 2 jF\ as determined from 
Eqs. (HU) and flH7|) of the text with tp 1 oc (M 2 + (0.6 GeV) 2 )" 3 - 5 and cp 2 fixed by 
Eq. fT8]) . The constituent mass is m = 0.3 GeV, and the bound-state mass is M = 
0.94 GeV. The choice of ip 2 corresponds to the Chung-Coester ansatz and nominally 
implies a constant asymptotic behavior for the ratio, but here is well fit by Q~ ' 25 . 

The arguments by Ralston et al. [H] and Kroll ji2], in favor of QF 2 /F\ being 
asymptotically constant, reduce to a discussion of the following ratio of wave function 
matrix elements: _ 

El = tyifo) (44) 
(See for example Eq. (5) of Ref. |H] and Eq. (11) of Ref. [12] •) They arrive at this 
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ratio by including an L z = 1 wave function in the contribution to F 2 , so that there 
is an overlap with the L z = wave function. However, the L z = 1 wave function 
^_i/ 2 also contributes to F\ in an overlap with itself, as given above in Eq. ()33)1 . This 
introduces an additional term of order Q 2 in the denominator and makes Q 2 F2/F\ 
asymptotically constant. 



5 The Wick— Cutkosky model and dependence 

The Wick-Cutkosky model is based on the ladder approximation to the Bethe- 
Salpeter (BS) equation in <p 2 x field theory. The LFWFs in this model for the 
J = L = and J = L = 1 bound states have been computed in jsJJjj. We shall 
show below that if the wave function is represented in four-dimensional form, the 
scalar components depend on Mq as expected on general grounds. Indeed, for J = 1 
the BS function reads |2S H3] 

*ix(q,p) = ^(?W?,p), (45) 



gi (z,M)dz 



[m 2 — M 2 /4 — q 2 — zp-q — ie) 4 ' 



where gi(z, M) is the spectral function satisfying a specific differential equation [2H 
SS]- We will not need the explicit form of g\. Here q = (k a — k b )/2, q = L _1 (-u)g has 
the sense of relative momentum in the rest frame given by p = 0, and L _1 ( , y) is the 
Lorentz boost with v = p/e p . In an arbitrary frame the spherical function in f)45|) can 
be replaced by 



M[|)--y£e? ) (p)^ (46) 



where eft'(p) is the spin-1 polarization vector. 

Using the relation (jSJ between the LFWF and the BS function and Eq. (JIEJ), we 
find the LFWF in the form 

^ix = -\^-el x \p)r, r = k^ + —^ (47) 

V 47T ^ UJ-p 



where 



2^(1 -2x,M) 
3(M 2 -M 2 fx 2 {\ -x) 2 ' 



9l (l-2x,M) g[(l-2x,M) 

^ 6(M 2 -M 2 ) 2 x 2 (l-x) 2 3(M 2 - M 2 ) 2 x(l - x) ' 1 } 

The prime in g[(z,M) means differentiation with respect to z. We see that the 
components (pi and (f2 depend on Mq, as expected. A similar result holds for the 
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L = state: 

* = g(l-2*,M) 

^ 2 v ^r(l-x)(M 2 -M 2 ) 2 ' 1 j 

where g(z, M) is the corresponding spectral function. It does not depend linearly on 
M . 

The representation (|4"7}l for the Wick-Cutkosky wave function is analogous to the 
representation (j2J) for the Yukawa model. These results for the LFWFs of the bound 
states of the Wick-Cutkosky model can also be obtained by computing the instant- 
form wave function and boosting to infinite momentum, as in Weinberg's P z — > oo 
method |%4] . 

Note that the LFWF (|4"7|) in the constituent rest frame has the form 



^u(M) = fi(k 2 ,n-k)kY* x ^-J + / 2 (^n.fc)Y 1 * A (n), (50) 

where the scalar components /i^ can be expressed in terms of tpi^ from (|47j) . This 
representation is analogous to the representation (|T3*|) for the Yukawa model. Since 
the BS function ()45|) describes a state with angular momentum J — 1, the correspond- 
ing LFWF (|50j) . which is derived from the BS function without any approximation, 
definitely has the same angular momentum. It is an eigenfunction of the angular 
momentum operator (JBJ) (omitting the spin operator 



6 Consequences of the CC ansatz for three-body 
light- front wave functions 

The simplified three-quark LFWF of the nucleons introduced by Coester and Chung, 
and used in [TH] for a model of the proton form factors, has the same spin structure 
as the quark-scalar-diquark model — two quarks form a spin-zero diquark, and the 
spin of nucleon is determined by the spin of the third quark. The wave function is 
symmetrized relative to the permutations of the quarks according to SU (6) flavor-spin 
symmetry. For example, in the constituent rest frame, the Pauli-spinor representation 
of the three-quark wave function corresponding to the state of zero spin and zero 
isospin of the quark pair 12 is 

^(12,3) = X (S=0) (1,2) X t 35 ^(/=0) (1)2) ^ra^ 

where x^ S=0 \^-i 2) = Xl^yX^ 1S ^ ne spin-zero wave function of the two quarks and 
£( 7=0 )(1 ; 2) = ^ Tl ia y SX 2 i s ^ ne isospin-zero wave function. The nucleon spin-isospin 
wave function is obtained by symmetrization: 

^(1, 2, 3) = ^(12, 3) + ^(23, 1) + ^(31, 2). 
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Using the Fierz identities, we may transform it to the form 



1p = —7E=1ps[3 + (v12-V3n)(t 12 -T 3N )}. (51) 

This form is totally symmetric with respect to spin and isospin. We have introduced 
here the symmetric momentum-dependent part ips- The factors in (joTj) should be 
understood as 

1 = (xljVyXlj (xl 3 X* N ), {viyVzn) = (xl^yXlj (xI/Xctn), 

and similarly for the isospin part. The wave function (|5T|) is just the ansatz used 
in [HI. 

In the covariant Yukawa model, the general form of the two-parton LFWF (fT3*|) 
contains two components, in contrast to the one-component CC ansatz ()12j) . The gen- 
eral form of the three-quark nucleon wave function contains sixteen components [3B] 
in contrast to the one-component ansatz (JSTj) . Thus the CC ansatz (Jo"Tj) is even a 
stronger constraint than in the two-body Yukawa model. 

As in the two-body case (jHJ), the three-body wave function ()51|) can be recast in 
the four- dimensional form [3*H] 

if) = —=c x c 2 c 3 c N {3[u(k 1 )A + j 5 U c u(k2)\ [u(k 3 )A + u N (p)\ 
V 72 

-[uik^A+Y^Ujlih)] [u(k 3 )A + -f^ 5 A + u N (p)}(T 1T T 3N )} , (52) 

where A_ = v 2 ~^° an d U c = 727o is the charge conjugation matrix. The coefficients 
are c\ = 1/ ' yjm + etc. The projection operator A + is again given by Eq. 0; how- 
ever, the two-body constituent momentum is replaced by the three-body momentum 

V = h + k 2 + k 3 , Ml = V 2 . 

In the constituent rest frame, where V = 0, we have Vq = Mo, and the wave function 
(jHg) reduces to (JST]) . 

The second term in (JSg) does not contribute to the proton form factor [33 UHj • 
The first term is factorized. We suppose that the photon interacts with the third 
quark. The first factor [u(ki)A + , -f5U c u(k2)] is the same for F\ and and, therefore, 
it does not change their ratio. Since only the second factor \u(k 3 )A + UN(p)] gives 
different contributions to F\ and F 2 , we rewrite the wave function in a factorized 
form 

ip ocips u(k 3 )A + u N (p). (53) 
This coincides with the CC-constrained Yukawa wave function (jHJ). 
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The form factor calculation using the three-quark wave function determined by 
(JoTj) and ()53|) differs from the two-body calculation in Sec. 0] only in the use of three- 
body kinematics. The result is 



m 2 ) 

F2{q 2 ) 



1 

2 

M 



x 3 (l - x 3 )(M 2 + M'l) + 2x 2 3 M M^ + 2x 3 m(M + M' Q ) - (1 - x 3 ) 2 Q 2 
xG(l,2,3,Q 2 )A (54) 



2(1 - x 3 )m + x 3 (l - x 3 )(M + Mq) - ^(M Q - M^) 2 (M + M' Q ] 
xG(l,2,3,Q 2 )D. 



(55) 



Here D is the three-body phase volume, and G(l, 2, 3, Q 2 ) is a function of the variables 
of all three quarks and of Q 2 . We do not need the explicit form of G. Mq is the 
effective three-body mass, which depends on R 3 ± — |(1 — x 3 )q±, whereas Mq depends 
on R 3 ± + ~(1 — x s )q±. Because of factorization of the wave function (|53j). these 
formulas (found by direct calculation) coincide, after evident changes of notation, 
with the corresponding expressions (j41J) and (}4~2*|) for the two-body form factors for 
the CC ansatz. 

At Q 2 = q 2 oo, both M' Q and Mq tend to Q. For the leading term we find 



F 1 oc ...Q 2 
F 2 oc ...Q, 



(56) 



which reproduce the ratio QFi/F 2 = const. As before, this asymptotic behavior is a 
consequence of the specific constraints on the LFWFs. 

On the other hand, we have mentioned in Sec. 0] that scalar and pseudoscalar 
couplings result in the asymptotic ratio F 2 /Fi ~ 1/ log(Q 2 /m 2 ) and vector coupling 
in F 2 /Fi ~ log(Q 2 /m 2 )/Q 2 . These couplings are not really related to a hypothesis 
of predominance of diquarks, but simply represent different spin structures of the 
nucleon wave function. The total number of these structures is sixteen j3H]. Some 
of these other structures, in addition to the vector one, may also contribute to the 
asymptotic behavior F 2 /Fi ~ \og(Q 2 / m 2 ) / Q 2 . 



7 Conclusions 

The explicitly Lorentz-invariant formulation of the front form provides a general 
method for determining the general structure of light-front wave functions. We have 
also used the fact that the angular momentum of a bound state can be defined covari- 
antly within the Bethe-Salpeter formalism: The LFWFs for an n-particle Fock state 
evaluated at equal light-front time a = u-x can be obtained by integrating the covari- 
ant Bethe-Salpeter functions over the corresponding relative light-front energies. The 
resulting LFWFs are eigenstates of the kinematic angular momentum operator ©• 
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The result is that the LFWFs are functions of the invariant mass squared of the con- 
stituents Mq = (YJ k^) 2 and the light-cone momentum fractions Xj = kfu/p-u, each 
multiplying spin-vector and polarization invariants involving u 11 , where uj = (ujq, uj) 
is the four-vector determining the orientation of the light-front plane uj-x = 0. 

We have presented the structure of LFWFs for two and three-particle bound states 
using the explicitly Lorentz- invariant formulation of the front form pQ. As examples 
we have given the explicit form of the LFWFs for spin-0 and spin-1 eigenstates of 
the nonperturbative eigensolutions of the Wick-Cutkosky model, as well as examples 
of spin-1/2 states constructed using perturbation theory. For example, the LFWF of 
a spin-1/2 system composed of spin-half and spin-zero constituents has the general 
form 



where the ifi are functions of the square of the invariant mass and the light-front 
momentum fractions. The orbital angular momentum prefactor in the constituent 
rest frame is proportional to uj x k-S. 

An important test of the LF computations is light-front invariance — although the 
LFWFs depend on the choice of the light-front quantization direction, all observables, 
such as matrix elements of local current operators, form factors, and cross sections, 
must be independent of We have computed the large momentum transfer behavior 
of the ratio of Pauli and Dirac form factors of the nucleon using the exact relation 
for spacelike current matrix elements in terms of LFWFs. The dependence of the 
invariant mass squared implies that hadron form factors computed from the overlap 
integrals of LFWFs are analytic functions of Q 2 in agreement with dispersion theory, 
the PQCD analysis of Belitsky, Ji, and Yuan j2Hj, and conformal arguments jH], as 
well as with the form factor ratios obtained using the nonperturbative solutions to 
the Wick-Cutkosky model. We have also shown that a fit to the Pauli to Dirac form 
factor ratio incorporating the predicted perturbative QCD 1/Q 2 and logQ 2 asymp- 
totic dependence describes the recent Jefferson laboratory polarization transfer data 
well. In contrast, we have shown that the LFWFs introduced by Chung and Coester 
to parameterize the static and low-momentum properties of the nucleons correspond 
to the spin-locking of a quark with the spin of its parent nucleon, together with a 
positive-energy projection constraint. These extra constraints lead to an anomalous 
linear dependence of the LFWFs on the invariant mass of the constituents and an 
anomalous dependence of form factors on Q rather than Q 2 . 

The CC construction of relativistic light-front wave functions was introduced to 
represent the properties of the nucleons in the low momentum transfer domain. How- 
ever, there are a number of difficulties with extending the CC form to the high 
momentum transfer domain: 

(1) If one applies the CC ansatz to a bound state of a spin-half quark and a scalar 
(the Yukawa model), the quark is constrained to have the same spin projection as 
the bound state nucleon, S* = S*, when one uses the conventional BD spinor repre- 
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sentation. Thus the only orbital angular momentum allowed by the CC constraint 
is the kinematical angular momentum arising from the lower components of the BD 
spinor. The spin-locked CC constraint does not allow for the full degrees of freedom 
of a relativistic system. 

(2) If one compares the CC ansatz to wave functions generated in perturbation 
theory, the net result is to replace the bound-state mass M in the numerator of 

the LFWFs by the invariant mass of the constituents M = \ E?=i Xi Z * ■ This 
replacement leads to the anomalous growth of the CC LFWF at large transverse 
momentum. For example, if one applies the CC ansatz to the QCD quark splitting 
function, the ultraviolet logQ 2 behavior of q(x, Q 2 ) will derive from two sources: the 
standard behavior in k± arising from perturbative QCD plus the presence of a term 
Mq oc k± due to the CC ansatz. The presence of the latter source would destroy 
DGLAP evolution in Q 2 . 

(3) The LFWFs of the J = and J = 1 bound states can be obtained ex- 
plicitly in the Wick-Cutkosky model. The form factor ratios of the spin-1 system 
obtained in this non-perturbative analysis is given by a quadratic Q 2 dependence and 
by logQ 2 |25j . The application of the CC ansatz leads to terms in the wave function 
which are linear in M rather than the quadratic dependence of the explicit solutions. 

(4) In the case of the simple Yukawa theory, the effective interaction in the CC 
model has the form HJ S = gip q Aip(f) where A is a non-local positive-energy projection 
operator. The presence of the projection operator A conflicts with the usual relations 
obtained from crossing and particle-antiparticle symmetry. For example, consider 
the electroproduction amplitude j*P ~ * Q + <?<?) where, for simplicity, the qq diquark 
can be taken as the Yukawa scalar cf). The Born amplitude with quark exchange in 
the t channel has the form 

M^ p _ q(f> (p-q,t,q 2 ) ex x i[>(x, k ± ) 

t m q 

where 

t-m\ = x(M 2 -M 2 ), 

2 

and x = XBj = The CC ansatz introduces a linear term in Mo in the electropro- 
duction amplitude. If we now use s — > t crossing to obtain the process 7*g — > pep, the 
presence of a linear term in Mq in the LFWF gives a contribution to the amplitude 
M^g^p^s, t, q 2 ) which is proportional to at fixed momentum transfer t. This 
anomalous Regge behavior corresponds to fermion exchange in the t channel, which, 
however, is not present in this amplitude. 

Again, we emphasize that these difficulties concern the extrapolation of the CC 
ansatz to the asymptotic region. These problems do not appear in the original 
work ^H] where the CC ansatz was only applied to static nucleon properties and 
to form factors at relatively small momentum transfer. 

Light-front wave functions are the fundamental amplitudes which relate hadrons 
to their fundamental quark and gluon degrees of freedom. We have shown how one 
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can exhibit the general analytic structure of light-front wave functions, including 
states with nonzero orbital angular momentum. A key element of this analysis is the 
use of the explicitly Lorentz-invariant formulation of the front form where the normal 
to the light-front is specified by a general null vector cj m . The resulting LFWFs are 
functions ip n {xi, k±i) of the invariant mass squared of the constituents Mq = (J2 k^ 1 ) 2 
and the light-cone momentum fractions Xi = kfUj/p-uj, which multiply invariant pref- 
actors constructed from the spin-matrices, polarization vectors, and in the case 
of nonzero orbital angular momentum. The LFWFs corresponding to definite to- 
tal angular momentum are eigenstates of a kinematic angular momentum operator 
and satisfy all Lorentz symmetries of the front form, including boost invariance. We 
have illustrated these properties using known nonperturbative eigensolutions of the 
Wick-Cutkosky model for nonzero angular momentum. The dependence of LFWFs 
on the invariant mass squared implies that current matrix elements and hadron form 
factors are analytic functions of Q 2 in agreement with dispersion theory and pertur- 
bative QCD. We have also shown that a model incorporating this analytic property 
and leading-twist perturbative QCD constraints is consistent with recent data for the 
ratio of proton Pauli and Dirac form factors determined by the polarization transfer 
method. 
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A caret indicates the inner product of a four-vector with the 7 matrices, so that 



Appendix 



We use the convention a 1 * 1 = a°±a 3 , which gives a-b 
use the 7 matrices in the Dirac representation 



|(a + 6 + a b + ) — a±-b±, and 




(57) 



k = fc/y". 

The solutions of u(ki)k 
spinors, are given by 



1 — 



mu(ki) and pu{p) 



Mu(p), in terms of the BD 



u 
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u a {p) = Je p + M\ g.p \x a , (5? 



e p + M 



where x a is a two- component spinor, Ek = yk 2 + m 2 , and e p = \Jp 2 + M 2 . In the 
constituent rest frame, where p + lot = k\ + k 2 = 0, we introduce the variables A;, n, 
ki = k, and u3 = nc^o- We also find 

_^M? — M 2 M 2 + M 2 

» = —cur = — n , e D = . 

y 2M p 2M 

The value of r here was obtained by squaring the equality p + our = k\ + k 2 , which 
gives r = (Mq — M 2 )/(2u-p), and by using u-p = u-(ki + k 2 ) = uj M . In this way, 
we find the BD spinors in the constituent rest frame to be written as 




a-k 



Ek + m 



M + M 



= vsn- ( M -M)*.ffK' (59) 

The factor in (171) is transformed as 

^ ~ Mo (7 ° ~ n ' l] ~W [n-a -1 J ' 

We also introduce 

uj-ki if n-k \ 1 n-k 

x ~^p~~2y ~7k~J~2~w ' 

Substituting these expressions into Eq. (jJJ, we reproduce the wave function (|T3j) with 
the components /i, f 2 given by Eqs. ([Hj) . 
The light-front spinor is defined as |10j 

which gives 



(A*) = X tai (k° 1 +m + a 3 (a-h), -fife? - m)a 3 - a-h) . (61) 

%k\ + fc?) V 7 



The BD and LF spinors are connected by the following unitary relations (at n\\z): 
( u + ^(p) \ 1 ( (p + + M) -p R \ ( u LF+ ^(p) \ 



u'Hp) J ^/2p+(p° + M) \ P (p + + M) j \ u LF -2(p) 

u LF+ Hp) \ _ 1 / (p+ + M) p R W u+3(p) 

u lf ~Hp) ) ~ ^/2 P +(pO + M) V -p l (P + + M ) ) V «"i(p) 



(62) 
(63) 
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